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Abstract 

The seasonal polar caps of Mars can be used to test the equivalence principle in 
general relativity. The north and south caps, which are composed of carbon dioxide, wax 
and wane with the seasons. If the ratio of the inertial to gravitational masses of the caps 
differs from the same ratio for the rest of Mars, then the equivalence principle fails, 
Newton’s third law fails, and the caps will pull Mars one way and then the other with a 
force aligned with the planet’s spin axis. This leads to a secular change in Mars’s along- 
track position in its orbit about the Sun, and to a secular change in the orbit’s semimajor 
axis. The caps are a poor Eotvos test of the equivalence principle, being 4 orders-of- 
magnitude weaker than laboratory tests and 7 orders-of-magnitude weaker than that 
found by lunar laser ranging; the reason is the small mass of the caps compared to Mars 
as a whole. The principal virtue of using Mars is that the caps contain carbon, an element 
not normally considered in such experiments. The Earth with its seasonal snow cover can 
also be used for a similar test. 
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I. INTRODUCTION 

Mars has an atmosphere composed primarily of carbon dioxide. Each Martian 
year a significant fraction of the C0 2 atmosphere freezes out in the form of polar caps [1, 
2]. These caps, which wax and wane with the seasons, can be used to test the equivalence 
principle in general relativity. 

The basic idea is the following. Assume for the moment Mars has only one polar 
cap. If the equivalence principle fails, then the gravitational mass of the cap attracts the 
rest of Mars with a gravitational force which differs in magnitude from Mars attracting 
the cap’s inertial mass. It follows that Newton’s third law fails and the planet will self- 
accelerate. The absence of a measurable self-acceleration indicates the equivalence 
principle holds, at least within the limits of error. The spirit of this astronomical test is 
thus similar to that of Bartlett and Van Buren [3], who used the heterogeneity of the lunar 
crust and the Moon’s lack of observable self-acceleration to make a stringent test of the 
equivalence principle. 

The qualitative details of the polar cap test are shown in the schematic diagram in 
Figure 1. Mars orbits around the Sun. For the purposes of illustration, Mars is shown to 
have an obliquity (axial tilt) of 90°. The Sun heats up Mars’s north pole when it is over 
that pole (bottom of the figure). The northern polar cap shrinks over time while the 
southern cap grows. This creates a net gain of CO? in the southern hemisphere (right) and 
a net self-acceleration (if there is one) along the spin axis, as shown by the thick arrow. 

The arrow opposes the motion, as shown in the figure, if the effect is such that the 
southern cap pulls on Mars more than Mars pulls on the cap. When the Sun shines on 
Mars’s south polar cap (top), that cap shrinks and the northern cap grows, and once again 
the self-acceleration opposes the motion (left) for the sign of the effect adopted in the 
figure. The end result is a net negative along-track acceleration when averaged over one 
revolution of Mars about the Sun; the semimajor axis of the orbit shrinks with time. 

The thick arrows are shown in Figure 1 as opposing the motion. However, 
because the sign of the self-acceleration is not known a priori without a theory, the 
arrows could just as easily point the other way. In this case, the orbit would expand 
secularly with time. The lag between the maximum insolation and the minimum size of 
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the cap is taken to be 90° in Figure 1 for ease of illustration; the lag for the real Mars is 
less than 90° (see below), but not zero. 


II. EQUATIONS OF MOTION 

The body of Mars will be taken to be a sphere. Each seasonal polar cap will be 
assumed to be a spherical cap centered on the pole, and with a uniform mass density 
which increases and decreases with time. The surface density of the northern polar cap 
will be On - win/An, where nnj is the mass of seasonal north polar cap, and An - 2k R 2 (1 
- cos d N ) is the cap’s area, where 0 N is the colatitude of the edge of the cap. Likewise, the 
southern seasonal cap has the corresponding quantities Os = ms/ As and 6s. In the 
following 6n and 6s are each taken to be constant. Thus a cap’s edge does not advance or 
retreat, as observed. The surface densities On and o s are constant in latitude or longitude, 
but vary with time. In other words, at any given time the CO 2 cover in each cap has a 
constant thickness, but the thickness changes with time. These assumptions will introduce 
an error whose magnitude is estimated below. 

The gravitational surface mass density ON g of the northern seasonal cap will exert 

a force 


df = 


GMp Ng dA 
R 2 


on the inertial mass Mj of Mars, where G is the universal constant of gravitation, R is 
Mars’s radius, and dA = R 2 sin 6 d8 dX is an element of area, with 6 being colatitude and 
A being longitude. The force df is directed radially outward from the center of Mars. The 
total component of force along Mars’s spin axis is then 

/ •h n r in „ GM . m v „ 

o f g sm0cos 6dAd0 = — ^ r - A -K(0 jV )p 
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where p is the unit vector in the positive spin axis direction, m Ng is the gravitational mass 
of the northern cap, and 


K{B n ) = 


1 - cos 2 d N 
4(1 -cos Q n ) 


( 1 ) 


The body of Mars will pull on the northern cap with a force Fn = GM g ms, K(6s)/R 2 ) in 
the direction opposite to fy, where ms, is the inertial mass of the northern cap, and M g is 
the gravitational mass of the rest of Mars. The acceleration of Mars will then be 
approximately 


a N P 


A-fA ram, 

\ r 


K{6 n ) p 


due to the northern cap alone, where 


r = 


m 


Ng 


m 


Ni 


M, 


( 2 ) 


Thus the test involves the CO 2 caps and the crust-mantle-core of Mars, whose 
compositions are quite different. The equivalence principle fails if F * 0. Note that the 
equivalence principle can still fail if ( ms g /ms ,) * 1 and ( M g /M t ) * 1, but ( ms g /mm ) = 
(M g /Mi); in such a case r - 0 and the failure will be invisible in this test. 

Analogously, the southern cap will give an acceleration as 


AsP “ ' 


fs-Fs 


M 


P = 


FGin 


p / 


R 2 


K(6 S ) p • 


in the direction opposite to on- The net acceleration of Mars along its spin axis will then 
be 
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«,vsP = K ~ a s )P = -~r\ m m K i e N) ~ m Si K{6 s ) Ip (3) 

A 

If the quantity in brackets is nonzero and F # 0, then Mars will self-accelerate. 

If all of the mass in a cap is concentrated at the pole, then K = 1 . Observations 
indicate that each polar cap extends to a maximum of about 35° in latitude from their 
respective poles [4], In this case = K(ds) = 0.91, so that (1) is fairly insensitive to 
cap size. In the following, it will be assumed that K{6 m) = K(6s) = K, so that (3) becomes 


a /vs P ^5 )P 


rGK 

R 2 


K - m si)P 


(4) 


and in the numerical calculation below K is taken to be 0.9. 

The northern and southern cap masses vary with time. The principal terms are 

m Np = c m sin {L - d m ) + c m sin {1L - 0 iV2 ) (5) 

m Sp = c S i sin (L - 6 S! ) + c S2 sin (21 - (5 52 ) (6) 

where L = co+f is the areocentric longitude of the Sun measured in a Mars-fixed frame, 
with m being the L of perigee and /being the true anomaly of Mars’s orbit about the Sun. 
Because Mars’s orbit is fairly eccentric with an eccentricity e = 0.09,/ does not increase 
uniformly with time. However, the mean anomaly M does increase almost uniformly with 
time. To average equations over time to get the secular effects on Mars’s orbit, it proves 
convenient to express L in terms of M. To first order in e [5] 

sin / a sin M (7) 

cos f= -e + cos M (8) 


sin 2/s -2e sin M+ sin 2 M 


(9) 
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cos 2/s -2e cos M + cos 2 M ( 1 0) 

Using these in the expression for mn P in (5) yields 

mm = -e cm (cos <5,vi sin oj - sin byi cos co) 

+ cm (cos dm cos co + sin <5 jvi sin ca) sin M 
+ cm (cos by | sin to - sin byi cos to) cos M 

- 2e Cm (cos 6m cos 2 to + sin dm sin 2 to) sin M 

- 2e cm (cos dm sin 2to - sin dm cos 2 to) cos M 
+ cm (cos dm cos 2 to + sin dm sin 2 to) sin 2 M 

+ cm (cos dm sin 2 to - sin dm cos 2 to) cos 2 M (11) 


with an analogous expression for ms,. 

The evolution of Mars’s orbital semimajor axis a with time t is given by 


da 

dt 


2 

nr( \-e 2 ) m 


[a( 1 - e 2 )S + (er sin f)U) 


( 12 ) 


where r-a( 1 - e 2 )/(l+ e cos f) is the Mars-Sun distance, n = Mars’s mean motion about 
the Sun, U is Mars’s acceleration directed radially outward from the Sun, and S is the 
acceleration in the orbital plane and is perpendicular to U, so that S is nearly the along- 
track acceleration for modest orbital eccentricities [6], Using (7)- (10) in (12) yields 


J ry 

— ~-(S + eScosM + eUsinM) (13) 

dt n 


to first order in e. 

Let r = cos L x + sin L y be the unit vector in the direction of U, where unit vector x 
points to the right in the plane of the orbit as shown in Fig. 1 , and y points to the top of 
the figure. The unit vector t = - sin L x + cos L y points in the direction of S. Unlike what 
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is shown in Fig. 1, Mars’s spin axis p = + sin 0y + cos 0 z lies out of the orbital plane, 
where 0 = 25° is the present obliquity of Mars and z is the unit vector normal to the 
orbital plane. Then S - (p *t) a ms and U= (p *t) ay.y, where 

(p «t) = sin 0 cos L sin 0(-e cos co- sin co sin M+ cos co cos M) (14) 

(p • r) = sin 0 sin L ® sin 0 (-e sin co + cos co sin M+ sin co cos M) . (15) 

by (7) - (10). By (1 1) and (14), and the acceleration 5 becomes 
5 « (G/Xsin ©HR 2 ) {-cm sin 8m + csi sin <5 91 

- 2 e[c m sin (co - 8m) - csi sin (co - fe)]} (16) 


ignoring the periodic terms and retaining only the secular terms. It turns out that by (14) 
and (15) eS cos M+ eU sin M- 0 to order e in (13), so that 


da ~ 25 
dt n 


(17) 


III. RESULTS AND DISCUSSION 

The following numerical values are used for Mars: R = 3.39 x 10 6 m, e = 0.09, n = 
1.06 x 10 ' 7 s' 1 , 0 = 25°, and L= 251° [7], while G = 6.67 x 10 ' 11 m 3 kg' s' 2 . Also, K is 
taken to be 0.9, as mentioned above. Smith et al. have recovered the coefficients and 
phase angles in (5) - ( 6 ) from the Mars Global Surveyor spacecraft orbiting Mars [2]; 
which are given in Table 1. These authors find that the cycles vary little from Martian 
year to Martian year. Using all these numbers in (16) and (17) yields 

S= 3.5 x 10 ' 9 Cm s ' 2 (18) 


and 
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da/dt = 2A x 10 8 Tracy' 1 . (19) 


One source of error in the above estimates is the value for K. This is because the 
carbon dioxide caps are assumed to have a fixed size, and only thicken or sublime. The 
thickness of the real caps vary both spatially and with time. However, K is fairly 
insensitive to size; even if 9n = 50°, K(9m) goes down to only 0.82. Hence it seems likely 
that setting K(9n) = K(9 S ) = K = 0.9 produces an error in the treatment of the caps of no 
more than -10%. 

Another source of error is the neglect of the Mars atmosphere. As the caps wax 
and wane the atmosphere thins and thickens. The atmosphere “feels” the topography, 
changing the C0 2 mass distribution and essentially creating another surface mass layer in 
addition to the caps. The atmospheric effect is about 10% that of the caps [8] and will be 
ignored here. 

It remains to estimate P Krasinsky and Brumberg find that the the solar system, 
as measured by the Astronomical Unit, may be expanding by 15 ± 4 m cy' 1 , for which 
they find no satisfactory explanation [9], Attributing all of the expansion to Mars’s polar 
caps in yields a value of | P | < 7.1 x 1 O' 8 by (19). 

The Earth’s changing snow cover will also generate a polar cap- type effect, which 
can complicate solar system tests, since the Earth is often used as one leg of a test. The 
Earth’s seasonal snow cover has a maximum mass of about 10 15 kg [10], which is about 
the same as the Martian CO? caps. The radius of the Earth is almost twice as great as 
Mars’s in (4), and the Earth’s mean motion about the Sun is also almost twice as great as 
Mars’s in (17). If the value for P for the Earth’s H?0 - (crust-mantle-core) system is 
comparable to the value of r for Mars’s C0 2 - (crust-mantle-core) system, then da/dt for 
Earth is -15% that of Mars, and a closer look at the Earth’s polar cap effect appears to be 
warranted. The Earth will not be as simple as Mars, because of the distribution of the 
snow on the continents, and the fact that sea level will go up and down as the snow cover 
changes. 
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Other solar system tests of the equivalence principle are also not considered (e.g., 
[11]). Given these complications and the absence of a more detailed study, it seems 
reasonable to increase the limit on Tby a factor of 1 0 to 

|r|s7.1xl0~ 7 (20) 

This limit is 4 orders-of-magnitude greater than that given by laboratory experiments, and 
7 orders-of-magnitude greater than that given by the Moon (see, e.g., [12]), so that Mars 
is a poor Eotvos experiment by comparison. 

Mars does have the virtue of testing the equivalence principle with a novel 
combination of elements. The caps are almost entirely CO 2 , so that the caps contain 
carbon. The composition of Mars’s body is not entirely certain, but a recent model has 
the crust and mantle consisting of Si0 2 (44.4% by weight), MgO (30.2%), FeO (17.9%), 
and AI2O3 (3.1%), plus 4.4% other elements; while the core consists of Fe (77.8%), S 
(14.2%), and Ni (8.0%) [13]. Carbon versus these other heavier elements is not usually 
used in tests of the equivalence principle. 
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TABLE I 


The coefficients and phase angles for the north and south cap masses. 


=(-1.534 ±0.028) x 10' 5 kg 
c N 2 = (-0.486 ± 0.029) x 10 15 kg 


5 m =43.64° ± 1.07° 
5nz = 47.88° ± 3.34° 


cs\ = (-3.058 ± 0.032) x 10 15 kg 
c S2 = (-0.917 ± 0.033) x 10 15 kg 

dsi =223.84° ±0.61° 

Ss2 - 52.85° ± 2.02° 
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Figure 1 

Schematic of the Mars equivalence principle test. Mars orbits the Sun. Mars’s spin axis is 
shown as lying in the orbital plane for the purposes of illustration. The polar caps (white) 
wax and wane with insolation. When one cap is larger than the other, there is a self- force 
(thick arrows) if the equivalence principle fails. The planet self-accelerates, leading to a 
secular change in the size of the orbit. For the case shown below the orbit shrinks. If the 
thick arrows pointed the other way, the orbit would expand. 




